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An original dispersion relation between the stationary coherent nonlinear optical responses by cur- 
rent and polarisation is obtained. The dispersion relation provides a new complimentary tool that 
can be employed to study light-induced charge transport models and facilitate experimental data 
analysis. It is shown that the origin of the coherent current and the dc-polarisation induced in a 
noncentrosymmetric crystal under illumination is related to the theory of the Berry phase and can 
be represented in terms of the renormalised geometric potentials. This renormalisation originates 
from the extra phase difference acquired by a carrier in the light field on the quantum transition 
between the electronic bands. The gauge invariance of the corresponding expressions for the current 
|' and the polarisation is demonstrated. 
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I. INTRODUCTION 



Investigations of the interaction of an incident light with the electronic subsystem of a noncentrosymmetric crystal 
have exposed a number of time-independent effects of the second order in the amplitude of the external field. One 
of these effects is the anomalous bulk photovoltaic effect (APE), which consists of the generation of a directional 
macroscopic current in the closed circuit of a spatially homogeneous semiconductor [H-IH . Another related phenomenon 
is the optical rectification effect (ORE), which implies that the crystal becomes polarised Q. Both of these effects 
are stationary second order responses to light perturbation on the summarised zero frequency - the APE along the 
current and the ORE along the polarisation. The mechanism of the APE consists of two channels of different origin: 
a ballistic channel Q and a coherent channel 043- The latter channel is independent of the carriers' scattering 
relaxation times and is determined by the interference of quantum transition amplitudes @, which is related to the 
^ i ■ effective shift of the carrier wave packet in real space [l(J • 

In the late 1980's, we ascertained for the first time a special relationship between the current and the dc-polarisation. 
This relationship is given by the duality transformation (lT| - [l3j 

CN i P{lu) = r^jM (1) 

i> . 

which, before passing to the limit r\ — > 0, allows us to express the ORE polarisation (P) through the coherent current 
(J) of the APE for a given frequency uj of the incident light. Here, the parameter 77 , which is a reciprocal averaged 
0^ ' relaxation time, mathematically plays the role of an intermediate variable that would usually be taken to zero in the 
final calculations in nonlinear optics. Equation ([1} has been obtained within the framework of the density matrix by 
considering the equation of motion for the coordinate operator. In theoretical practice, the expression has allowed us 
to calculate the ORE polarisation by using the result of the summation over the intermediate electronic states 

In general, equation (J]) reflects the relationship between dissipative and non-dissipative processes in the system. 
The product rjP(u>) can be interpreted as the unique source of coherent APE current; in the limit r\ — > 0, the 
product creates selection rules for electron- hole generation induced by the absorption of photons with optical transition 
energies. It thereby provides steady-state charge transport in a noncentrosymmetric semiconductor. Because the 
processes of energy exchange between the light wave and the crystalline medium depend also on the frequency of the 
electromagnetic field, one can expect that the duality transformation in Eq. (p} corresponds to a wider dispersion 
relation that connects spectral distributions of the APE and ORE. In the present work, we derive such a generalised 
relation. As a practical tool, the relation will be especially useful for experimentalists because it provides a simple 
procedure for relating the spectral representations of the relevant polarisation and current to one another. 

On the other hand, an interconnection between the coherent APE current and geometrical features [3] of the Bloch 
states [Hldll has not been completely realised. Moreover, the key issue of the gauge invariance of the current remains 
unexplained. 
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The present study is a direct extension of our previous work [l7j and puts an emphasis on clarifying both the 
macroscopic and quantum-mechanical connections between the APE and the ORE. There are at least two main 
objectives for this work: 

(i) We wish to relate the current and the polarisation from the perspective of dispersion correspondence in addition 
to Eq. ((T|). Such new correspondence, established with respect tow, is of considerable importance because it would 
allow direct expression, within the framework of experimentally measurable frequency dependencies, of the coherent 
part of the APE current through the cfc-polarisation induced by the light wave. This link, in turn, makes it possible 
to experimentally represent one effect in terms of the other. 

(ii) In light of the application of Berry phase [l4[ effects in various fields of physics [l8[ , we wish to obtain a detailed 
understanding of the geometrical origin of the coherent APE current and the ORE polarisation. This is a natural 
step in analysing these macroscopic physical quantities because, from a quantum mechanical point of view, they are 
determined by interference processes in the space of Bloch states. By making this connection, we recall that the 
general conception of the crystalline polarisation is based on the Berry phase jm , [l9| . It is also evident that when 
seeking insight into the geometrical features it will be important to establish the relevant gauge invariance. 

II. THE DISPERSION RELATION BETWEEN THE COHERENT APE CURRENT AND THE ORE 

POLARISATION 

Second order interaction processes between light and matter have been investigated in a number of our works @, 
Ell Ell [13; the results from these works will be used in our present theoretical analysis. In the context of the duality 
of Eq. (TTJ), the functional formulation of the coherent APE current can be written as a sum of Lorentzian functions: 

J ^) = ^ E E J ^E / d 3 kC nl (k,e^ V )^(k) (2) 

e— ±1 (3= x,y n,l 

where 

C ^^) = (E n (k)-E l (l) + eu J r + V i> (3) 
*£f (k) = |p^(k)| 2 [/„(k) - /,(k)] Re{F^(k)} , (4) 

FZ?(k) = [flgck) - K n (k)] + \ d ln ^ (k)) . (5) 

Here Ip is the intensity of an elliptically polarised light wave, C is a distinct constant, p^Ak) with n^l and i?" n (k) are 
the matrix elements of the momentum (p = — i d/dr) and coordinate (i d/dk) operators on the periodic Bloch ampli- 
tudes Mnk(r), the indices a and /3 denote the Cartesian components, and f n (k)=f(E n (k)) are the Fermi occupation 
factors of the electronic states with band energies E n (k). The response by the coherent APE current is represented 
by = lim^o J a (oj) ; correspondingly, the ORE polarisation is P(u>) = Iim^o V Ja{u) ■ 
Observe that from the relation 
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where V denotes the principal value, one can readily write 



V 

=p iTrS(x) 



P(u) = - lim / dx !L P(x) . (6) 

7r »;^o J _ oc [x — ujy + rj z 

By considering the ORE polarisation and the coherent APE current as a dual pair with correspondence given by ([TJ , 
one finds 

P(") = T^v3{x) (7) 
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where integration is defined in the sense of Cauchy-Hadamard (e.g., [2CJ, |21|). Alternatively, there is another repre- 
sentation of (0 which is related to a Cauchy principal-value integral by 



I d f 



+ °° dx . . 

— — j(x). (8) 



If one defines the Hilbcrt transform of a real function of the real variable j(oj) as (e.g., |22i. [23j) 

1 f + °° dx , x 
H{j(u)} = -V j(x), 9 

7T J _ OQ X-UJ 

it is seen from Eq. (j8]) that the ORE polarisation is a derivative with respect to u> of the Hilbert transform of the 
coherent APE current. By using the property dW{j (uj)} / duj — H{dj (u) / duj} , Eq. (jHJ can be rewritten as 

PH = i/ M Af). (10) 

7T J _ OQ X — LO dx 

With the help of Eqs. ([9]) and (fTU| the ORE polarisation can be regarded as the convolution of the frequency derivative 
of the coherent APE current with a power function of the form l/(x — u>). 

The above results (Eqs. , (JSJ) , and (1101) ) constitute our desired dispersion relation, which connects the frequency 
dependencies of the coherent APE current and the ORE polarisation. Therefore, once we know a spectral represen- 
tation of one quantity the above equations give the relevant representation of the other. In fact, this provides a new 
complementary tool for the study of steady-state nonlinear optical effects. A more practical motivation for these 
studies lies in their technological importance. For example, the phenomenon of optical rectification is of considerable 
interest because it is related to the ability to generate terahertz radiation (e.g., [24( ) . 

Because P(oj) = P(—oj) , relation (0 can be also written as 

P(o J ) = -i dx * + " j(x) = - j-W dx^^j(x). (11) 



77 J -oo ( x ~ uj2 ) 2 77 du 



— oc 



As a simple illustration of Eq. (fTT|) . the spectral frequency profile of P(ui) for a ^-shaped coherent APE current 
= jo [S{uj - luq) + 5(uj + tdo)] is given (for lu ^ uj ) by 



III. THE GEOMETRIC PROPERTIES 

Let &f (\s) be the phase of the interband matrix element of the momentum operator: p>L(k) = |p^,(k)| e 1 ®'™^ . 
The first derivative of the amplitude d\pf n (k)\/dk a does not contribute to the integral in Eq. @ and, correspondingly, 
Eq. ([5]) can be written as 

Ftf(k) = [iZg(k) - K n (k)} + -g^- . (13) 

The first term in the square brackets of Eq. (|T3|) describes the difference in the positions of the electron charge for 
the initial and final states. By using the definition of a gauge potential (a Berry connection) in the k-space for the 
band n as A„= R„ n , one has 

Ff„ = A n ~ A l+ grad6f„ . (14) 

As seen from Eqs. - (J5J) and (fl"4")) . the expression for the current in the geometrical sense is formulated in terms 
of vector potentials and not fields, and so the relevant gauge invariance must be verified. Clearly, the current is an 
observable physical quantity and should be independent of the freedom to change the phase of the Bloch wave function 
arbitrarily at each point in k-space. 

The gauge symmetry of the current is simply stated from Eq. ([5|) by considering the quantity F^ n to be invariant 
under a local gauge transformation: u n k(r) — > e l Vn ^> M nk (r) where the quantity <£>„(k) is an arbitrary phase of 
the Bloch amplitude. According to the sense of the Berry potential, we can call such a gauge geometric. As a 
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direct consequence, we reveal the physical role played by the gradient term in Eqs. ([5]) and (fl4|) as follows: the 
term maintains the geometric gauge invariance of the coherent APE current and the ORE polarisation, and hence 
eliminates the arbitrary non-gauge contributions caused by phases of the Bloch amplitudes. In view of Eq. (fT4"|) . this 
allows us to conceptually treat the phases of the matrix elements of the momentum operator as compensating fields. 
Now that we have described the roles of both components in Eq (fT4|) , it is possible to gain better insight into the 



basic properties of the quantity Ff n . First, by defining sf n as half of ©f n , we can rewrite Eq (fl4| as 

Ff„ = A r „ - Af (15) 

where A^ = A„ + gradsf n and A[ = A; + gradS^, , respectively. Furthermore, we call attention to the point that, 
due to the Stokes theorem, the calculation of the difference between the Berry phases $„ and $/ , which are acquired 
by the wave functions of the Bloch initial ip n u(^) and final V'/k( r ) states, does not involve the gradient terms 



*j = jd\-[A n - A,] = idl-Fl = Idl-[A r n - A[] . (16) 



A key observation that directly follows from Eq. (fT6|) is that there is a simple way to rewrite the quantities A^ by 
employing the definition of the Berry connection 

K = << k (r)Md/dk|< k (r)) (17) 

where the matrix elements are calculated on the Bloch amplitudes renormalised by multiplying by the phase factor: 
U 



(r) = u„ k (r)exp(-iSf n 



Such a representation of the phase-renormaliscd Bloch states, in which the relevant gradient terms move into the 
phases of the Bloch amplitudes, leads to a new understanding of the role played by the phases of the interband matrix 
elements of the momentum operator and, therefore, allows us to draw several conclusions: 

(i) From Eqs. (p~5|) and (ITTl) . one sees the geometrical meaning of Ff n as the quantity representing a difference of the 
renormalised potentials of the n and I electronic bands. 

(ii) The phase factor exp(— iSf n ) provides a direct argument that an extra nontriviality appears in the geometry of 
the Bloch states of the perturbed system, a nontriviality that is of geometric phase origin (a similar extra phase, often 
called "statistical", takes place in studies of the statistics of quasiparticles entering the quantum Hall effect [25j). 
Geometrically, the allocated phase factor is what incorporates the effect of the other band in the Bloch amplitude of 
the given band through the half-phase of the momentum matrix element between these bands. This, in turn, implies 
that it is also possible to consider the quantity sf n in the sense of phenomena such as the Aharonov-Bohm effect [IB] 
as the phase difference adiabatically acquired by the electron on the transition between the n and / bands. 

(iii) These phase differences cannot be ignored because they are employed as compensating fields to ensure that 
equations for the coherent APE current and the ORE polarisation are independent of the choice of the Bloch amplitude 
phase. Thus, the geometric gauge invariance of both physical quantities is provided by a nontrivial superposition of 
the difference between the "bare" Berry connections and the interband phase difference acquired in the electron 
transitions. 

We must emphasise that the coherent photoinduced current of the APE and the efc-polarisation of the ORE are 
completely defined by the gauge invariant contributions of the renormalised Berry connections in the space of the 
Bloch states' quasimomenta. Because they are determined by the geometry of this space and not by the particular 
dynamics of the Bloch electrons, neither effect naturally depends on the specific scattering mechanisms of the carriers. 
For the present case, the characteristic geometric features are associated with the interference of quantum transitions 
between the intermediate electronic states. Therefore, using the terminology of the anomalous Hall effect [27], such 
contributions can be called intrinsic contributions. 

In this situation, the meaning of the relaxation parameter ?y, which, like the Berry phase, is strictly related to the 
adiabatic evolution of the system, can also be deepened. One may say that the parameter limits the time to collect 
the necessary Berry phase change during the light period. 



IV. DISCUSSION 



The final question on which our study can shed light concerns the interrelation of the analytical and geometrical 
properties of the coherent APE current and the ORE polarisation. Following a causal formulation of the coherent 
nonlinear photoinduced effects [l7| , consider the evolution of the electronic subsystem of a crystal in the external light 
field. According to the concept of time ordering in the dynamic quantum multi-electron system [28 30] and by using 
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the results of Ref. [17| , an expression giving the stationary second-order response for an arbitrary operator D can be 
rewritten in terms of the Heaviside step functions 8(i) as follows: 



SM<£»= r 2 e~^ J db x J dt 2 e^ 0(t - ti)0(*i - t 2 )([[D(t), V(h)], V(h)])\ o ■ (18) 

Causal-like sequencing of interactions (time ordering) in the time evolution of the system is imposed here in an 
explicit form through the products of the Heaviside step functions. In Eq. (TT51) . V(t) is the perturbation caused by 
illumination, the tilde means the interaction representation and index "0" selects the stationary part of 5^ <D> . 

A valuable property of a representation like Eq. (fT8|) is the possibility to employ a time sharing scheme 0, [3(3| in 
which the time ordering is separated into two terms - the uncorrelated (T unc ) and correlated (T corr ) parts - by using 
the following representation of the function 9(ti — £2) : 

9(ti - ta) = - ( ^ + sign(ti - t 2 ) ) . (19) 

Clearly, the two contributions tend to complement each other. The first contribution, T unc , corresponds to an 
independent time approximation where the V(tj) interactions are disconnected among themselves in time. The second 
contribution, T corr , connects the V(tj) at different times [2{|. By employing the Fourier transform of <d(ti — £2), the 
given scheme (|19[) can be supplemented with two parallelisms [28l . |30j 

T unc — > lim — I 2 = tt 8(E - E ± oj) , (20) 

tj— >-q [Eq — E ± ujy + r] z 

Tcorr ► lim r-= rr = — (21) 

77^0 (Eq -E±Lj) 2 +7] 2 Eq — E ±Ll) V ' 

where E is an excited state in the electronic subsystem. 

In order to see whether the coherent APE charge transport and the ORE polarisation can be categorised according 
to their time-ordering properties, we shall consider separately each contribution prescribed by Eqs. (fT9| - (|2"Tj) . As 
a first case, we discuss the situation driven by T unc (i.e., the time- independent part of 0(£)) in which the sequence 
of the interactions V(tj) is unimportant. As mentioned in Ref. [2{|, the second order process can be represented in 
this approximation by two independent one-step processes, which correspond in Eq. (|20p to the delta function. At 
the same time, as seen from Eq. ([3]) at 77 — > , the one-photon interband absorption processes are responsible for the 
generation of free carriers. Because they are dependent on the relevant amplitudes, these processes maintain no phase 
information inside the area occupied by electrons. This implies that we need to invoke the system's other degrees of 
freedom in order to discover the asymmetry responsible for the rise of the coherent APE current In quantum 

mechanical systems, one available degree of freedom belongs to the special geometric properties in the k-space of the 
Bloch states [HI]. As we have seen above, the microscopic origin of a macroscopic coherent current is related to an 
interplay of the Berry geometric potentials and the interband phase differences. In a medium without a centre of 
inversion, the sum of these quantities provides gauge-invariant polar distributions in individual electron transitions 
between the Bloch states in momentum space. 

Therefore, one can conclude that the light-induced stationary coherent APE charge transport, being of geometric 
origin, does not need any special time correlations between electrons. This result is not surprising because one should 
expect that time-correlating processes cannot be involved in the case under consideration due to the absence of the 
time-reversal invariance associated with the bulk macroscopic current. Moreover, from the quantum mechanical point 
of view correlation in time, which is initially an attribute of the incident photons due to Bose statistics, is destroyed 
along with the photons themselves, which are absorbed in order to generate the electron transitions. 

A diametrically opposite situation takes place in the second case when the time-dependent part of <d(t) in Eq. (fT9]) 
is considered. This case includes the time-correlating processes in the time evolution of the electron subsystem excited 
by the external electromagnetic field. In a mathematical context, the case is given by the principal value contribution 
in Eq. (|2"Tj) and driven by T corr . From Eqs. @, ©, and (f2"Tj) . one sees that the light-induced ORE polarisation 
P(w) can be associated with the effect of the time correlation term, T corr . This result implies that the interactions 
of a light wave with the electron subsystem at various times must interfere to polarise a nonlinear medium. The 
need for such interference can be described in a more pictorial way. The external perturbations affecting the electron 
subsystem lead to local changes of the electron density distribution that, in turn, create imbalances with respect to 
the corresponding density of the positive background. In the case of the time correlations, these imbalances will be 
accumulated in a sequential and dependent manner in order to involve a sufficiently large number of valence electrons 
in virtual transitions. The induced polarisation is, thus, a result of a new macroscopic redistribution of the electronic 
charge. 
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Following from the above discussion, time ordering is the result of a nonlinear effect due to the specific correlations 
caused by the interactions of a light wave with the electronic subsystem. This leads to interesting consequences with 
respect to the two physical quantities related by the duality transformation in Eq. ([lj : one turns out to be connected 
with the time correlation while the other is not (at least in the second order approximation). Simultaneously, both 
of them are caused by the geometric features in the k-space of the Bloch states. In addition to the explanation given 
in the introduction, this observation allows us to treat the factor 77 as a driving parameter, which shifts the time 
correlated and uncorrelated contributions amongst themselves at r\ — > 0. Next, analogously to the above treatment 
of fTJ), the physical meaning of the dispersion relation in Eq. (JT]) can also be completed. One can say that this relation 
interconnects the time correlated contributions with the uncorrelated contributions with respect to u. Obviously, such 
an add-on to the explanation of (J7J should have greater generality because, for example, a similar treatment exists 
for a much different situation (i.e., fast collisions of atoms with matter and light) [29j . 

In summary, we focused our attention on the analytical properties and the geometric meaning of the second- 
order photoinduced effects in noncentrosymmetric crystals. We established an original dispersion relation, which 
connects the spectral representations of the coherent APE current and the ORE polarisation. It was shown that 
both phenomena are encoded by the renormalised Berry connections (potentials). The renormalisation of the Berry 
connections is determined by the interband phase differences in the quantum transitions between the Bloch states and 
provides the underlying gauge symmetry in the k-space. We compared properties of the coherent APE current and 
the ORE polarisation with respect to a physical picture in which the electron subsystem is considered to be correlated 
in two different ways. It was determined that one way takes into account the time-correlating processes created by the 
light wave while the other way includes the correlations related to intrinsic geometrical features of the Bloch states. 
As a result, we observed that the ORE polarisation appears to be the combined result of both ways. Conversely, the 
coherent APE current is solely associated with the second way. 
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